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We summarize some recent progress in the understanding of the statistical properties of 
cosmological billiards. 
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1. Introduction 

A remarkable achievement of theoretical cosmology has been the construction, by 
Belinski, Khalatnikov and Lifshitz (BKL), of a general solution to the 4-dimensional 
v acuu m Einstein equations in the vicinity of a spacelike ( "cosmological" ) singularity 
12121^ HI They found that this solution exhibits a never-ending oscillatory behavior, 
with strong chaotic properties. They could describe in detail the statistical prop- 
erties of this never-ending oscillatory behavior by approximating the Einstein field 
equations (near the singularity) by a system of ODE's for three variables a, &, c 
('anisotropic scale factors'), namely 
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where dr = —dt/{abc), and by approximately reducing the continuous dynamics 
of a, 6, c to a sequence of discrete maps. The crucial discrete map introduced by 
BKL relates the 'Kasner exponents' pa, Pb, Pc describing the (approximately linear) 
r-evolutions of the three scale factors a, 6, c during two successive 'epochs' (i.e. two 
successive segments of the dynamics ([1]) during which the influence of the right-hand 
side is negligible) . Later studies have refined the description of the statistical prop- 
erties of the chaotic BKL oscillations, notably by introducing and studying more 
complete discrete iteration maps (involving several real variables), and notably by 
a two-dimensional discrete map 

Separately from the work of BKL, and, with a different aim and motivation, Misner 
realized that generic Bianchi IX homogeneous cosmological rnodels have a "very 
complex singularity" ^ a reformulation of this dynamics led to the sim- 

ple picture of a point moving on a Lobachevsky plane and reflecting upon fixed 
billiard-type cushions. This led Chitre to remark that the dynamics of the system 
point is ergodic and mixing, with unique invariant Liouville measure. The problem 
of relating the statistical properties of the discrete BKL map to the invariance of 
the Liouville measure, in the continuous billiard dynamics, a la Misner-Chitre has 
been considered in some detail by Kirillov and Montanim^. 

The description of cosmological singularities in term of billiards in (higher dimen- 
sional) Lobachevsky (or Lorentzian) spaces has recently received a new impetus 
from the discovery that the billiard chambers corresponding to many interesting 
physical theories can be identified with the "Weyl chambers" of certain (infinite- 
dimensional) Lorentzian Kac-Moody algebras ^I^, This has raised the con- 
jecture that, hidden below the BKL "chaos", there lies a remarkable "Kac-Mo ody 
symmetry" , akin to the duality symmetries of supergravity and string theories 
HZl, More precisely, in Lorentzian space, the Lagrangian C for these models reads 



L = iGabP'^p'^ - V{P) (2) 

where Gabd^^dfi^ = J2a i^P^ f - (^„ dp^f, and where the potential term V{P) = 
^ Cyie"^™'**^'^-' determines the billiard walls wa{P)- In the case of pure grav- 
ity, the only walls that can appear are either the gravitational walls w^abciP) ~ 
- - + P ot: the symmetry walls w'J^iP) = - In 4 = 3 -I- 1 
dimensions, the usual BKL billiard table is obtained by considering the three grav- 
itational waUs wl2z = 2/3^, a23i = 2/?^, 0312 = 2/3^. This 'big bilhard' corresponds 
to considering a diagonal metric^ described by the three scale factors a = , 
b ~ e^^ , c = . The more general case of a non-diagonal metric corresponds to 
a 'small billiard' enclosed by the gravitational wall 0123 = 2/3^, together with the 
two leading symmetry walls = 0^ ~ and ^23™ — P^ — P'^- Projecting the 

3-dimensional Lorentzian (/?) space on a Lobachevsky plane allows one to obtain 
some useful representations of the billiard dynamics. 
We focus here on a few results from Ref. [181 



January 13, 2013 16:50 WSPC/INSTRUCTION FILE proceeding 



About the Statistical Properties of Cosmological Billiards 3 

2. Descriptions of billiards 

The dynamics on a billiard table consists of free-flight evolution and of bounces on 
the boundaries of the billiard. The trajectories describing the free- flight evolution are 
called (Kasner) epochs, while a set of trajectories taking place in the same corner 
of the billiard table is called a (Kasner) erg0. Six kinds of eras can be defined, 
according to the corner where an era takes place, and the orientation of the first 
epoch of the era. All the information necessary to describe the dynamics can be 
encoded in maps that take into account collisions on given walls. 




Fig. 1. The billiard tabic in the unit-disc rep- 
resentation. 



u - 1 H M 

Fig. 2. The billiard table in the Poincare 
half-plane representation. 



There are two complementary descriptions of billiards: the unit disk and the 
Poincare half plane, sketched in Fig. [l]and Fig. [21 respectively. 
The unit disk is the 2-dimensional version of the general n-dimensional 'ball' confor- 
mal representation, geometrically realized [within the (n-l-l)-dimensional Lorentzian 
/?-space] by stereographically projecting (from the 'South Pole' 75, i.e. a center of 
projection located on the past unit hyperboloid) the future unit hyperboloid onto 
a (n-dimensional) hyperplane passing through the origin in /3-space. This represen- 
tation is powerful in outlining the symmetry group of order six that characterizes 
the 3 -|- 1-dimensional version of this model. 

The Poincare model can be obtained by a suitable geometric inversion of the ball 
model. The Poincare half plane offers a less symmetric representation, but it offers 
a convenient parametrization both of the dynamics and of the oriented endpoints of 
the geodesies. Each collision can be described by a linear fractional transformation 
of the variables (of the plane) : 

^, ^_az + 6^ ad-6c=l. (3) 



^consistently with the definition given in Eq. (5.4) of Ref. |4] 
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For each given wall, there exists an application of the form ^ which transforms each 
geodesic incident on the considered wall in the reflected geodesic. When parametriz- 
ing a geodesic by the two real numbers parametrizing its end points on 
the 'absolute' (see Fig. [2]), the transformation ([3]) also applies, separately, to w_ 
and u+. Note, in particular, that the single real parameter m+ parametrizes the 
Kasner exponents of the corresponding Kasner epoch. This m+ parametrization of 
the Kasner exponents actually coincides (as emphasized by Kirillov and Montani) 
with the usual BKL parametrization, namely 



u 



Pi'-i^) = - ..2 , . , 1 . (4a) 



2.1. Integral invariants for billiards 

Several integral invariants can be constructed for billiards. Their role is to provide 
conserved quantities (that can be used in the statistical analysis) charactering the 
dynamics. 

The usually considered energy-shell Liouville measure is a three-form obtained by 
reducing to the hypersurface H — E the four- form ojpQ A Wp^, where ijJplj = 
dq^ A dpi — dt A dH{q,p, t) is the (invariant) Poincare-Cartan 2-form. However, one 
can also construct another invariant measure by directly restricting the Poincare- 
Cartan two-form on any transverse section of the Hamiltonian flow. In the case of 
the Poincare plane (asymptotic) Hamiltonian for billiards, 

H^{u,v,Pu,Pv) = {pI +pI) + Voc{i{u,v)), (5) 
the transversally reduced Poincare-Cartan two-form reads 

^reduced ^ ^ ^ ^ dp.],=,o„,t = A dp^ = 2^^^^. (6) 

The two- form ^ is then found (both by general arguments, and by explicit calcula- 
tion) to be invariant under the discrete 'billiard map', which connects one collision 
to the next one. Each billiard map is given by Eq. ([3]) and is therefore of the form 

/ au-\-+h 
U+ ~ ^t:j. 

The invariant two-form ([6]) is linked to the Liouville 3-form by VL)^ = i^^edMced ^ '^^^ 
where ds is a phase-space parameter measuring the length along each geodesic seg- 
ment. 

3. Big billiard hopscotch game 

We describe the dynamics of the hig billiard as a hopscotch game at different levels, 
such that different properties can be stressed upon. 
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The epoch-hopscotch description of the bilhard dynamics is obtained by dividing 
the complete bilhard table into six regions B^y, each denoting the portion of the 
reduced phase space where all the epochs going from the wall x to the wall y take 
place. Fig. [3] illustrates the 'jumping around' of successive epochs in the {u~,u~^) 
plane: from the epoch number 1 (from b to a) to the epoch number 5 (from b 
to c). Each transformation between successive epochs is given by linear fractional 
transformations of (u~,m+). 



Bba 










1 


.5 Bbc 


-1 


Bac J 
















Bcb 





Fig. 3. The epoch-hopscotch court. The six different colors (shades of gray) denote the six regions 
of the bilhard table, where the different kinds of epochs take place. 

We can also use the symmetry group S3 (of order 6) to quotient the dynamics 
of the big billiard. One way would be to consider a kaleidoscopic version of the big 
billiard dynamics in which the single "moving ball" of the billiard is augmented by 
its 5 (generically distinct) images under S3. This leads to a billiard game where 6 
(symmetry related) balls simultaneously move within the same billiard table, and 
(simultaneously) bounce on its bounding walls. A second way to look at the quo- 
tiented dynamics is to replace the latter kaleidoscopic phase-space point 
{q^-^^, 9(2): •■•'9(6);^'^^-P^■^^ ■•■'-Pi''''} by its unique representative, say g^ep, within a 
fundamental domain of S3, together with its corresponding momenta Pmp- When 
passing from the continuous billiard dynamics to the discrete billiard map from an 
epoch to the next epoch, the quotienting of the big billiard leads to a quotiented 
version of the epoch hopscotch game of Fig. [3l In the alternative, fundamental- 
domain, version of the quotiented dynamics we could replace each S3 orbit in the 
(w~, u+) plane by its unique representative located within, say, the box Bba- This is 
achieved by means of a suitable Kasner transformation. [Indeed, the six boxes B^y 
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of Fig. [3] are exchanged under S'3.] In that view, the discrete quotiented big bilhard 
would become a map from the box Bjja onto itself. If we ignore the u~ coordinate, 
we see that the law giving the successive values of the coordinate coincides with 
the law found long ago by BKL. This shows that the BKL discrete dynamics of the 
variable u is obtained by: (i) quotienting our more complete hopscotch dynamics 
by the permutation group S'3, and (ii) ignoring the W coordinate and identifying 
the BKL variable u with . 

The discrete epoch hopscotch map (both in its unquotiented and in its quotiented 
versions) leaves invariant the measure w = /^reduced ^ ^ addition, by 
marginalizing uj over u~ , this leads to the following invariant one-dimensional mea- 
sure for the usual BKL map on the u parameter (in the u > convention): 

wta{u^)du+ = . (7) 

As far as we know, this result has not been explicitly discussed before in the liter- 
ature. 

However, one cannot avail oneself of usual results of ergodic theory because the 
integral of w on its total domain (e.g Bba for the quotiented map Tia) is infinite 
(similarly the integral of Wbadu~^ is logarithmically divergent). In order to be able 
to use the usual results of ergodic theory, one needs a (discrete) invariant map pre- 
serving a finite measure. This can be achieved by considering the era hopscotch 
map (either in the unquotiented or quotiented forms). An era hopscotch dynamics 
is obtained simply by ignoring the intermediate epochs and focusing on the map 
transforming the (possibly quotiented) first epoch of an era {up,Up) into the first 
epoch {up , it^ ) of the next era. The unquotiented era-hopscotch court is sketched 
in Fig. m , and the explicit expression of the unquotiented era map T is listed in 
Table[TJ In the quotiented case, the resulting discrete (two-dimensional) map, say T, 
coincides with the map studied by Chernoff-Barrow and by Lifshitz-Khalatnikov- 
Sinai-Khanin-Shchur. When using a representative of the S3 orbit within the Bba 
box, the quotiented era map T is a map of Fba onto itself, i.e. the rectangular 
domain ~2 < Up < —1, < Up < +00. The explicit expression of the map T is 

T«l = + ± \ - 1- (8) 

Up - [uj.\ 

This maps leaves invariant the restriction of the two- form uj to the domain Fba] 
by contrast with the original measure uj on the full hopscotch court, this restricted 
measure has now a finite integral, namely J^^ ujp = 2 In 2. 

3.1. Recovering information about the unquotiented dynamics 

Quotienting the dynamics of the big billiard implies a loss of information with 
respect to the full unquotiented dynamics. We analyse here what kind of information 
is lost, and show how to recover this information. 
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Fig. 4. The era-hopscotch court. The six different colors (shades of gray) denote the six regions 
of the biUiard table, where the first epoch of the different kind of eras can take place. The asterisks 
indicate the one-epoch eras that form the simplest periodic configuration in the big billiard. 

Within the notation we have developed, the continued-fraction expansiontEl of the 
first u^) values of a (quotiented) era starting in the box Fi,a can be written as 

u^ = nba-l + [n2,n3,n4, ...] = [uba - 1; ?i2, ?^3, ?^4, •■•], (9a) 
Up^ -1- [mi, 1712,1713,...]. (9b) 

In terms of these decompositions pertaining to the first era, we can write the first 
u~,u~^ values of the iV-th era (with N = 2,3,...), i.e. the {N — l)-th iteration of 
the T map 



-N-l 



lit = ?lAr - 1 + [77N+l,nN+2,nN+3, 



-T-N-l - 

I Ur 



= -1 - [njv, n^v-i, . 



,n2,nba,mi,m2, ■ 



(10a) 
(10b) 



Furthermore, if we start an unquotiented era hopscotch dynamics in some specific 
region F^y, we can first map it to the reference region F^a by some specific 'Kasner 
transformation' k^y (listed in Table [5]) to get its 6a-representative, u^^°'\ namely 



+ [ba\ 



1 + [712,713, 



(11) 



''The notation [x] for a: G R denotes the usual integer part of x when x > (e.g. [it] = 3), 
and — [— < when a; < (so that [— vr] = —3); we define the fractional part of a negative 
number accordingly, e.g. {— tt} = —0.14.... Furthermore, we have denoted the first integer of the 



decomposition of uj, as n^^ 
starting with ui,. 



so that = [v 



1 denotes the length of the era 
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We see that the quotiented era dynamics encodes the lengths of successive eras, but 
ignores the fuller information contained in the original, unquotiented era hopscotch 
dynamics, namely the precise corner, {xjv,2/jv} (with xjv,t/Ar G {a, 6, c}) and orien- 
tation, {xn, Un), (i-e. xn — >• Un), of the first epoch of the A^-th hopscotch era, such 
that one could express the explicit result of iterating (iV— 1) times the unquotiented 
T map in terms of the simpler action of the quotiented T map, namely 



The missing information is recovered by considering the quantities 

p^ = e'^", with 9N^-D[kF]^iei+eie2 + - + eie2-eN), (13a) 
Vn = eie2-eN- (13b) 

where ej = (— )"^^^, rij denotes the number of epochs contained in the j-ih era, 
and D[kp] is the determinant of the first Kasner transformation k^y involved in 
()12|1 . The procedure determining k~^ is as follows; (i) starting from kxy and the 
continued-fraction decomposition of k^yU^^ = ni — 1 + [712,71.3, ...], one determines 
the ei, and the pjv and rj^] (ii) then pjv determines the rotation (in the unit-disc 
representation) between xy and the final corner xmtVn, and r/^r determines whether 
this corner is "parallel" {rj^ = 1) or "antiparallel" {rj^ = — 1) to the initial xy; (iii) 
finally, knowing the ordered corner [x^ ,y^), Table [2] determines the transformation 
k^N yN that maps it onto the ba corner. 



3.1.1. Periodic orbits 

Periodic orbits offer another way to contrast the unquotiented dynamics with its 
quotiented version. First, any m-periodic orbit of the one-dimensional BKL map 
Tbkl, parametrized by the special values of that admit a (regular) periodic 
continued-fraction expansion of the type 

u+ + 1 = [rii;7i2,7i3,--.n-m,ni,"2,"3,---] , (14) 

gives rise to a unique corresponding m-periodic orbit of the two-dimensional T map. 
Second, any m-periodic orbit of Tbkl and T can be lifted to some periodic orbit of 
period mp of the full, unquotiented billiard, according to the condition 

T"'Piu-,u+) = kP{u-,u+) = {u-,u+) , (15) 

where fc* is a suitable Kasner transformation (which depends on m and on the 
considered periodic orbit of T), and p is the order of k^. The order of a particular 
group element, such as fc,, is the smallest integer p such that fc» = fco! as a trans- 
position is of order 2, and a cyclic permutation, (123) or (321), of order 3, we see 
that the order p of fc* must be equal to p = 1, 2 or 3. In other words, the period 
in the unquotiented billiard can be up to three times larger than the one in the 
quotiented one. For example, the simplest 1-periodic orbit of T given by the golden 
ratio 7i+ + 1 [1; 1, 1, 1, ...] = (1 + \/5)/2 corresponds to a 3-periodic orbit of T, as 
shown in Fig. [5l 
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Fig. 5. The periodic succession of one-epoch eras of the big billiard in the unit-disc representation. 
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Table 1. The era-hopscotch transformations, which de- 
fines the big-biUiard map T- 
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Table 2. The Kasner tranfor- 
mations. 
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